clc;
clear;

load('EDC.mat');

% 提取涡旋数组
eddyArray = EddyDatasetClass.eddyArrayAddByTraj([], EDC.traj, -1);

% 取测试数组
eddys = eddyArray(10000:11000);

EAMSC = EddyArrayMatchSSTAClass(eddys);

% 求涡旋平均三维结构和平均半径
[SSTA, R] = EAMSC.mean_eddyMatchedSSTA_Radis();


% 求单个涡旋三维结构
eddy = eddys(1);
[eddyMatchedSSTA, Radis] = EddyArrayMatchSSTAClass.eddyMatchedSSTA_measure(eddy);

pou = squeeze(SSTA(31,31,:));
%% 注意需要更改82行  armorfilePath = 'H:\ftp\3D sst\月度平均armor滤波后数据\' 这个路径
